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Abstract 

Particle diffusion in a curved tube embedded in R3 is considered. We find the diffusion coefficient 
depends on tube's curvature. Diffusion coefficient is obtained in e (radius of tube) expansion. 
Physical interpretation of curvature dependent diffusion coefficient is given. 
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I. INTRODUCTION 



In the previous paper on surface diffusion with thickness embedded in R 3 established the 
curvature dependence of surface diffusion [1] . This theory might be allied to the molecular 



a. 



diffusion in lipid bilayer or applied to the reaction diffusion system |2[ 

However, the physical meaning of curvature dependence was not clear. To make clear 
this point we consider an one dimension less system, that is, the diffusion in tube embedded 



sional extension is given in {5]. (This quantum mechanical problem was originally discussed 
by many authors Q to solve the ordering problem of quantum mechanics on curved 
space.) In the case of thin tube (string) in R 3 , the dimensional difference is two and there 
appears additional degree of freedom to rotate around the central axis directions, so called 
torsion that makes calculation bit complicated. However the "local equilibrium condition" 
(explained later) makes the effect of torsion disappear, and we have the curvature dependent 
diffusion equation. The physical reason of its curvature dependence is discussed precisely in 
last section. 

II. FRENET-SERET EQUATIONS 

Let us consider a line specified by 

x(s), (I) 
where the parameter s is the length of line. The unit tangent vector is defined by 

ei = (2) 
as 

Another (normal) unit vector e* 2 is defined by equation. 

dei(s) _ 
ds = Ke 2 , I e 2 |= 1. (3) 

k is called the curvature of this line in R3. The reason is the following. We can set a 
circle which tangents to this curved line at point x(s). From figure 1 we can easily find two 
relations. 

e*i(s + ds) — e\(s) = e*2 d9, (4) 
ds = R d6, (5) 



in R3. The quantum mechanical version to this problem is given in [4|, and higher dimen 



FIG. 1. Curvature 




FIG. 2. Torsion 

where R is the radius of this circle. Then we find the equation ([3]) with identifying k = 1/R. 
Next we introduce another independent unit vector e 3 by 

e 3 = ei x e 2 . (6) 

Then we can identify the derivation of e*2 generally such as, 

— = a ei + r e 3 , (7) 
where a and r are unknown some functions. But a can be calculated as follows. 

_ de 2 _ dei 
a = ei • — = -e 2 • — = -«. (8) 
as as 



So we can replace (J7J) as 



— = -k e 1 + r e 3 . (9) 



ds 

Function r has a geometrical meaning. 
From the equation 

r = — ■ e 3 , (10) 

we see r means the rotation of normal vector e*2 around el (figure 2). So this function is 
called torsion. 
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We further take a derivative of e 3 . 



^ = /3e 1+7 e 2 , (11) 
with unknown function and 7. Both of them can be obtained as follows. 

n de 3 _ dei _ _ _ 

p = — ■ e x = — -7- ■ e 3 = -ne 2 ■ e 3 = 0. (12) 
as as 

and 



de 3 _ de 2 _ f . 

7 = ^- e2 = -^- e3 = - T - (13) 

So we obtain 



d ^ 3 - flA\ 

- = -r e 2 . (14) 

We call these three equations (j3J), (EJ), and ( fT4l as Frenet-Seret equations, and we have 
two geometrical quantities, curvature k and torsion r. 



III. METRIC IN TUBE 

We re-identify the one dimensional diffusion on curved line as the limitation process from 
three dimensional diffusion. We set the curved tube with radius e in three dimensional 
Euclidean space R3. Our particles can move only in this tube. We look for the form of 
diffusion equation in the limit of small e. The coordinates we use hereafter is the followings. 
(See FIG.3) 

X is the Cartesian coordinate in R 3 . x is the Cartesian coordinate which specifies only 
the points on center line. q % is the normal coordinate in a direction of e*j. (Small Latin 
indices i,j, k, ■ ■ ■ runs from 2 to 3.) s is the length parameter along the center line. Further 
by using the normal unit vector e*j, we can identify any points in tube by s,q 2 ,q 3 by the 
following equation j4]j5]. 

X(s,q 2 ,q 3 )=x(s) + q l e l (s), (15) 



where < | q | < e with | q | = a/ (q 2 ) 2 + (<? 3 ) : 



tei Tube 



Line / 








x// 




ee 3 


r = 



/' x(s) 

X{s, q 2 ,q 3 ) = x(s) + q 2 e 2 {s) + q 3 e 3 (s) 

FIG. 3. Coordinate 

From this relation we can obtain the curvilinear coordinate system in tube (c R3) by the 
coordinate g M = (s, g 2 , g 3 ), and metric G^ u . (Hereafter Greek indices /z, z/, • • • runs from 1 to 
3 with s = g 1 .) 



_ dX dX 
^ ~ dq» ' dq»' 



Each part of G^ u is the following. 



(16) 



Gu = 1 - W + (k 2 + r 2 )(g 2 ) 2 + r 2 (g 3 ) 2 , 

G12 = -rg 3 , 

G13 = rg 2 , 

G 23 = 0, 

G22 = C33 = 1. 



(17) 
(18) 
(19) 
(20) 
(21) 



Compared to the previous paper, we have non zero off-diagonal elements. This comes 
from the fact that dimensional difference between outer space R3 and inner space R\ is two. 
Now the following relations follow. 



2\2 



G = det(G llu ) = (l-Kq 2 ) 



(22) 



The inverse metric is given as, 



(1 - nq 2 y 



x 



-rq* 
-T 2 q 2 q 3 

,2\2 



\ 



( 1 rg 3 

rg 3 (1 - Kq 2 ) 2 + (rg 3 ) 2 
^ -rg 2 -r 2 g 2 g 3 (1 - /tg 2 ) 2 + (rg 2 ) 2 J 
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IV. DIFFUSION FIELD IN TUBE 



Let us denote 3 dimensional diffusion field as (f)^ 3 \ and Laplacian as A^ 3 ). Then we have 
the equation with normalization condition 

90(3) 



dt 



DA( 3 V (3) , (23) 
1 = J <P i3 \q\q 2 ,q 3 )VGd 3 q} (24) 

where D is the diffusion constant, and G = det(Cr^). Our aim is to construct the effective 
one dimensional diffusion equation from 3D equation above. 



^ = DA^^\ (25) 
1 = f (f,W( s )ds, (26) 



where is the one dimensional diffusion field, and A^ 6 -^ is unknown effective ID diffusion 
operator which might not be equal to simple ID Laplace Beltrami operator d 2 /ds 2 . 
From two normalization conditions, we obtain 

1 = / ^\ q \q 2 ,q 3 )VGd\ 

[ J dq 2 dq\<P {3) VG)} ds, 

(1) (s) ds. 




Therefore we obtain the relation, 

0«(a) = J cf>®VGdq 2 dq 3 . (27) 
We multiply \[G to equation fl23|) and integrate by g 2 , g 3 , then we obtain 

^1 = D J (VGA®)<f>®dq 2 dq 3 . (28) 
From the form of Laplace Beltrami operator 



A (3) = C~ 1/2 C 1/2 C^ 



our diffusion equation has form 



D I ^ 1/2G ^ (3) )^ 3 - ^ 



dt 

Here we suppose the "local equilibrium condition" such as, 

-^T = 0, i = 2,3. (30) 

Because in direction of e*2, e*3, equilibrium would be holds in a short time dt ~ e 2 /D. 
When our observation is given in time scale t satisfying t » St, we can always assume 
local equilibrium condition ([3"0l) . 

Then we have instead, 

= D^-(J G^G^dq'dq 5 )^ 

+ D{J ^(G^G*)dqW}§-/ 3) . (31) 

From (I2T|) and (I3"0l) we also have 

0« 



0( 3 ) = iV = /" v^gW = vre 2 . (32) 



So we obtain 



' G^G^dq'dq 3 )-^-^ 



dt 7re 2 ds J ds 



The explicit calculation gives the integral. 

J ^-(G^G^dq^q 3 = 0, (34) 



and 



f G^dqW = f^L 
J J 1-KQ 



2^3 
2 



Kq" 

- 2 {l + (f) 2 + 0(e 4 )}. (35) 



7 




Then we come to the result. 



FIG. 4. Curved point 



with the definition of effective diffusion coefficient, 



(36) 



D eff = D(l + (-y + 0(e*)). 



(37) 



The static solution has the form, 



0« = d + C 2 f 
Jo 



ds' 



Deff(s') 



(38) 



V. PHYSICAL INTERPRETATION 



We find the relation that, 



D, 



where 



< 



eff 



>- 



D < 



1 



1 — KCf 



>, 



vre 2 / 



dq 2 dq 3 



(39) 



Then we find the simple interpretation. Let us consider the point P on tube where the 
curvature is k. We chose two sections near P, and discuss about the length connecting these 
two sections. See FIG. 4. At the coordinate q 2 , the length between two section is given by 
As' = (1 — nq 2 )As. Therefore the curvature dependence of effective diffusion coefficient is 
given by the mean value of rate of length i.e., < As/ As' >. 

Next we consider the physical meaning. The diffusion coefficient is proportional to the 
mobility by Einstein's relation. The mobility is also proportional to the conductivity with 



D = a 

z — ' S ; 



FIG. 5. Flux of Tubes 



<72 



lid 



I' 



FIG. 6. another cross section 



unit length and unit cross section. In the following we calculate the rate of conductivity 
between bent and straight for the same length and same cross section flux. 

The conductivity of flux consist of N infinitesimal thin tubes is proportional to the sum 
of each conductivity of tubes. Each conductivity proportionals to the cross section Act, and 
inversely proportionals to the length Sj. See FIG. 5. So we have 



When the flux is straight with length As , each length is the same As = si = s 2 = • • • , 
and so we obtain 



(We consider the diffusion only in s direction since we used the local equilibrium condi- 
tion.) So we obtain the rate 





where a 



In the case of FIG. 4, we have 





Then we obtain the result. 
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Dbent — < Z q > D stra - 

1 - nq 2 

On the other hand as we have shown as R = case in the previous paper [l( , the curved 
surface with one direction has no curvature and other has finite curvature has anomalous 
diffusion coefficient such as, 

D eff = (f + eV/12)D. 

This relation also can be explained in this physical interpretation. In this case we consider 
the cross section as in FIG. 6. 

i i r /2 ( , , r w . 1 i 

-e/2 



e/2 

dg 2 {l + «;g 2 + (/€g 2 ) 2 + ---} 



•■> 



-6/ 

,2 ,2 



1 + eV/12 + • ■ ■ • (40) 
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